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Abstract 

We analyze in detail the discrete-time quantum walk on the line by separating the 
quantum evolution equation into Markovian and interference terms. As a result of 
this separation, it is possible to show analytically that the quadratic increase in 
the variance of the quantum walker's position with time is a direct consequence of 
the coherence of the quantum evolution. If the evolution is decoherent, as in the 
classical case, the variance is shown to increase linearly with time, as expected. 
Furthermore we show that this system has an evolution operator analogous to that 
of a resonant quantum kicked rotor. As this rotator may be described through a 
quantum computational algorithm, one may employ this algorithm to describe the 
time evolution of the quantum walker. 
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1 Introduction 



The study of computational devices based upon quantum mechanics, i.e. 
Quantum Computation, has drawn the attention of researchers in the last 
few decades [1,2]. The recent advances in technology that allow to construct 
and preserve almost perfectly quantum states, have opened the possibility of 
building useful quantum computing devices. However, relatively few quantum 
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algorithms that outperform classical ones have been found [3,4]. The classi- 
cal random walk is an example of stochastic motion that has found classical 
applications in many fields. The quantum version of this problem has sev- 
eral features which are markedly different from the classical walk [5,6]. As 
some classical algorithms are based on random walks, it seems natural to ask 
whether quantum random walks might be a useful tool for quantum compu- 
tation [7]. 

The discrete-time quantum walk on the line was introduced as a generalization 
of the classical random walk to the quantum world. Here, we will focus on the 
discrete-time quantum random walk on the line from a new perspective which 
emphasizes the role of coherence as the physical reason behind the striking 
differences found in the quantum version of the walk. We first briefly introduce 
the basic notions and notation relative to the discrete-time quantum walk on 
the line. Consider a particle that can move freely over a series of interconnected 
sites. The discrete quantum walk on the line is implemented by introducing 
an additional degree of freedom, the chirality, which can take two values: 
"left" or "right", \L) or \R), respectively. This is the quantum analog of the 
classical decision of the random walker. At every time step, a rotation (or, 
more generally, a unitary transformation) of the chirality takes place and the 
particle moves according to its final chirality state. The global Hilbert space of 
the system is the tensor product H s <g> H c where the Hilbert space associated 
to the motion on the line is H s and the chirality Hilbert space is H c . 

If one is only interested in the properties of the probability distribution, it has 
been claimed [5,8,9] that it suffices to consider unitary transformations which 
can be expressed in terms of a single real angular parameter 9. Let us call the 
operators that translate the walker one site to the left (right) on the line in 
H s as T_ (T + ), and \L)(L\ and \R){R\ the chirality projector operators in H c . 
We consider transformations of the form, 

U{9) = {T_ <g> \L)(L\+T + ®\R)(R\}o{I®K(9)}, (1) 

where K{9) = a z e l6ay , I is the identity operator in H s , and a y and a z are 
Pauli matrices acting in H c . The unitary operator U (9) evolves the state by 
one time step, 

Mt+l))=U(9)m)}- (2) 

One of the most remarkable characteristics of the quantum walk on the line is 
that it spreads over the line faster than its classical counterpart. In this work, 
we apply a general approach which leads to an physical insight of why this is so. 
To do this, we rewrite the evolution equation as the sum of two separate terms, 
one responsible for the classical-like diffusion and the other for the quantum 
coherence [10]. As we shall see, the terms responsible for the diffusion obey a 
master equation, as is typical of Markovian processes, while the other includes 
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the interference terms needed to preserve the unitary character of the quantum 
evolution. In Section 2, we review the decomposition of the evolution in these 
two terms. In Section 3, we obtain analytical expressions for the first and 
second moments of the probability distribution for these terms, in the case of 
the quantum random walker. Finally, in Section 4, we discuss our conclusions. 



2 Derivation of the Master Equation from the unitary evolution 



In a recent work [10], we have shown in detail how a unitary quantum mechan- 
ical evolution can be separated into Markovian and interference terms. This 
approach provides a new intuitive framework which proves useful for analyzing 
the behavior of quantum systems in which decoherence plays a central role. 
It is particularly suited to describe the evolution of quantum systems which 
have classically diffusive counterparts. 

The unitary evolution associated to U(9) can be decomposed into a Markovian 
term and an interference term. We begin by expressing the wave vector, as 
the spinor 

!*(*)>= E (M5V>> ( 3 ) 



where we have associated the upper (lower) component to the left (right) chi- 
rality and the states \m) are eigenstates of the position operator corresponding 
to the site m on the line. The unitary evolution for \ty(t)), corresponding to 
eq.(2), can then be written as the map 

a,i(t + 1) = a i+ i(t) cos 9 + b i+ i(t) sin 9 

bi(t+ 1) =a^i{t) sm9-b i - 1 {t) cos9. (4) 

Note that for the particular case 9 = 7r/4, the usual Hadamard walk on the line 
is obtained. The cases 9 = and 9 = it/2 are trivial motions not considered 
here. We define the left-distribution of position as Pmhit) = |om(^)| 2 and the 
right-distribution as P m R{t) = \b m (t)\ 2 . Then, the probability distribution 
for the position is P m (t) = Pmh{t) + Pmn(t) and these distributions satisfy the 
map 

P hL (t + 1) = P i+1 , L (t) cos 2 9 + P i+1 , R (t) sin 2 9 + i+1 (t) sin 29 

P itR (t + l) = Pi-iAt) sin 2 ^ + P,_ liR (t) cos 2 £-A-iW sin 29. (5) 

where the interference term in eq.(5) has been renamed as f3i = K [aib*] with 
^R,(z) indicating the real part of z. 
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It will prove useful to write equations (5) in the form 

oo 

Pis(t + 1) = E T ^siPji(t) + Pis(t) sin 29 (6) 

j=— oo 
l=L,R 

where s, I take the two values of chirality L, R, f3 iL = f3 i+ i and /3m = —f3i-\. 
The transition probabilities are defined as 



Tu + i t LL — Tu_i ji R — cos 2 9 
Tu + i t LR = Th-\,rl = sin 2 9 

Tijjs = otherwise. (7) 

Note that these transition probabilities satisfy the necessary requirements 
Tij tS i > and J2i, s Tij, s i = J2j,iTij tS i = 1. Now it is clear that if the in- 
terference term (3i S (t) in eq.(6) can be neglected, the time evolution of the 
occupation probability is described by a Markovian process in which the tran- 
sition probability (i,s) — > (J, I) in a time At = 1, is given by T^^i. As is 
characteristic of Markovian processes, the new position and chirality depend 
only on the previous values for position and chirality. Since the chirality is an 
auxiliary dimension introduced to implement the quantum walk, we focus on 
the evolution of the position distribution for the particle, Pi = Pil + P%r , 
which can be obtained from eqs.(5) as, 

Pi(t + 1) = [P t+1 (t) + P t _ 1 (t)]cos 2 9-P l (t-l)cos29 

+ [A+i(*)-A-iW]sin2#. (8) 

We have used the fact that P m (t — 1) = \a m ^i(t)\ 2 + \b m+ i(t)\ 2 , a relation 
that is a consequence of the map (4). Note that if the interference terms are 
neglected in eq.(8), the resulting evolution is Markovian in two time steps. In 
the general case, eq.(8) does not possess a continuum equivalent, even in the 
absence of interference terms. However, there are two special cases where it 
is possible to obtain continuum limitsfll]. One is if we take the increments in 
space (Arc) and the increments in time (At) in such a way that the velocity 
v — ^| and the ratio are kept constant when At — > 0; in this case one 
obtains the Telegraphist's equation [12]. 



dt 



Q2p Q2p 



dx 2 dt 2 



(9) 



In the second case, v 2 At is kept constant as At — > and 9 can take any value. 
Then the more familiar classical diffusion equation is obtained. 



dP d 2 P , , 

w = D ^ < 10 > 
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where, in the previous two equations, D is the diffusion coefficient, 

Note that this expression reduces to the classical value, D = 1/2, for the 
Hadamard walk. The inclusion of the interference terms in the continuum 
equations (9,10) is done by adding 2cot$|^ to their right hand sides. 



3 Moments for the position distribution 



The evolution of the variance, a 2 = M 2 — Mf, of the distribution Pi(t) is 
a distinctive feature of the quantum walk. It is known [14] that it increases 
quadratically in time in the quantum case, but only linearly in the classical 
case. We obtain the evolution of the variance analytically from the evolution of 

oo oo 

the first and second moments, defined as M 1 = J2 i P% an d M 2 = J2 i Pi 

i=—oo i=— oo 

respectively. The evolution equation for these moments is, from eq.(8), written 
as 

Mi(t + 1) = 2 cos 2 ^ Mi (t) - cos 29 M 1 (t - 1) -2 sin 29 ^A(f) 

i 

M 2 (t + 1) = 2cos 2 0[l + M 2 (t)] -cos20M 2 (t-l)-4sin20£;i/3 i (*). (12) 

i 

In differential terms, these equations become, 

d 2 M 1 n 9 n dM 1 A n ^ n . . 

+ 2tan 2 — - + 4tan0 = 



dt 2 dt 
t¥l + 2 tan 2 9 ^ + 8 tan £ ift(t) = 2 (13) 

(JjL (XL 



3.1 Decoherent evolution 



If the sums in eqs.(13) corresponding to the interference terms can be ne- 
glected, the evolution of the moments of the distribution is given by 

d 2 M x n 2n dM l n 
— + 2tan 2 — - = 
dt 2 dt 

d 2 M 2 n o n dM 2 n 

+ 2tan », = 2. (14) 
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The general solution for these equations is of the form 



Mi(t) = C u + C 12 e' 2ttan 8 

M 2 (t) = C 22 + -\- a t + C 2l e~ 2tt ^ e (15) 
tarn 9 

where Cu, C± 2 , C 2 \ and C 22 are constants. For times larger than r = cot 2 9, 
the transient exponential is negligible and the variance increases linearly in 
time with a slope given by eq.(ll), 

a 2 ^D(9)t. (16) 

This is consistent with the expected result for a classical random walk on a line. 
We emphasize that we have obtained this result neglecting the interference 
terms in the exact quantum evolution, but this result is valid for any 9. 



3.2 Unitary evolution 



If the interference terms in eqs.(13) are not neglected, the process is not Marko- 
vian but unitary. The solution for arbitrary 9 is cumbersome and not particu- 
larly illuminating. Therefore, in this subsection, we particularize our solution 
for the case of the Hadamard walk setting 9 = 7r/4 in the map (4). Using 
Fourier analysis, the solutions for the amplitudes ai{t) and bi(t) can be ob- 
tained [5]. For the particular initial conditions %(0) = 8jo and bj(0) = these 
solutions are 



aj(t) = 

hit) = 



i + (-i) j+t r* dk 



2 J- n 2tt 

1 ■■- {-iy +t /- dk 

-7T 2"7T 



1 + 



cos k 



vT+ cos 2 k 



D -i(w k t+kj) 



VT + COS 2 k 



-i(w k t+kj) 



(17) 



where sin Wb 



Jk = ^ rom these amplitudes, an expression for the variance 
can be obtained in the long time limit. The sums in eqs.(13) are evaluated, 
using the Fourier expansion of the delta function, with the result 



E m = a 

j=~oo 
oo 

£ jPj(t) = -At + B. 



(18) 
(19) 



j=-oo 



The numerical constants are given by A = [2 - V2j /4 and B = 1 - 5y/2/8. 
Even though we have obtained these results in the long time limit, for finite 
times they correctly describe the average trend. We have checked this fact 
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numerically, by computing the relative differences A A/ A = J2j Pj ~ A /A 

and AB'/B' = [EjjPj + At - b] /(-At + B) as a function of time. The 
result is shown in Fig. 1. 



Using eq.(12), the evolution for the first and the second moments can be 
expressed as the maps 

M^t + 1) = Mi(t) -2A 

M 2 (t + l) = M 2 (t)+4At+(l-AB). (20) 
The solutions for these maps are 

Mi(t) = -2 At + C 

M 2 (t) = 2At 2 + (1 - AB - 2A)t + C , (21) 

with C and C arbitrary constants. The variance a 2 = M 2 - Ml is then 

a 2 ^2A{l-2A)t 2 . (22) 

This result is consistent with previous those of previous works. In ref. [14], it 
is obtained through numerical simulations, in [5] it is found through Fourier 
analysis, while in [13], the same expressions are obtained through a summa- 
tion over different paths. However, it is important to emphasize that, from 
eqs.(18) and (22), a direct relation between a coherent evolution and the long 
time quadratic increase of the variance has been established. In a decoherent 
evolution A and B are negligible and the increase of the variance becomes lin- 
ear in time, as seen in eq.(21). This is a particular instance of the Markovian 
process discussed in the previous subsection. 



3.3 Generalized quantum walk 



In this subsection we present an alternative analytical approach which can 
be applied to the generalized quantum walk, i.e. for arbitrary values of 9. 
We start by rearranging the original map, eqs.(4), to uncouple both chirality 
components. The resulting independent evolution equations are 

Oi(t + 2) -Oi(t) =cos9 [a i+1 (t + 1) - a i - 1 (t + 1)] 

bi(t + 2)-bi(t) = cos 9 [b i+1 (t + 1) - b^it + 1)} . (23) 

For long times the left hand sides of the above equations can be approximated 
by time derivatives. The result can be expressed as 

2^ cos (24) 
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Fig. 1. The fractional difference between the asymptotic value and the finite time 
value of the expressions in eqs.(18) and (19). The quantities AA and AB' are defined 
in the text. 

where 6 stands for either chirality component, a» or fej. We define the effective 
time i! = —tcosO, then eq.(24) becomes 



2^ = 6-1(0-^(0, 



(25) 



which is the recursion relation satisfied by the Bessel functions. Thus, we can 
write the general solution as 



+ 0O 



&(t) = E(- i r^(o)^-K* cos#), 



(26) 



where 6(0) are the initial conditions to be used in the differential equation 
(25). These initial conditions are not necessarily the same as those to be used 
in the discrete map (23), because the approximation of a difference by a con- 
tinuous derivative does not hold for small times. The solution (26) with appro- 
priate initial conditions provides a good long-time description of the dynamics 
of the discrete map (23). In this context, long-time implies many applications 
of the discrete map. Note that eq.(26) provides the additional information [6] 
that the long-time propagation speed of the probability distribution is given 
by cos 9. 

It is worth mentioning that the long-time solutions for the spinor amplitudes, 
eq.(26), have the same form as the time-evolution of the amplitudes of the 
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Fig. 2. Fractional difference between the variance as obtained from the approximate 
expressions (29) and (30), and the exact value a 2 from the original Hadamard map, 
eq.(4). The quantity Ac 2 is defined in the text. The initial conditions used in eqs.(29) 
and (30) are (f°) = 0.70206 ({), (~^ 2 ) = -0.05963 (J) and zero otherwise. For the 

calculation of a 2 , the particle is initially at the origin with chirality (£°) = (J). 

kicked rotor in a principal resonance [15,16]. In fact, if K is the dimensionless 
strength parameter of the kicked rotor as denned in [17], the time evolution 
for these amplitudes is given by eq.(26) after the substitution 



cos 9 



K 
4np 



p integer. 



(27) 



We can obtain analytically the increase in the variance implied by eq.(26). 
The position probability distribution can be expressed as 

Pitt) = E(- 1 )" (,+, ' ) [ S *(0K(0) + &,(0)&|',(0)] Ji-i(tcos9) J^tcosO). (28) 

LI' 



The first and second moments of this distribution are 

Mi(i) = -t cos 9 E K [ a »(°) a i-i(°) + ^(0)^-i(°) 



+ Mi(0) (29) 



and 

M 2 (t) 



|i + E» [a,(o)ar_ 2 (o) + 6,(o)&?_ 2 (o)] } 

- t cos9 £(2Z - 1)3? fe(0)a;_ 1 (0) + 6i(0)6| , _ 1 (0)l + M 2 (0). (30) 
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respectively, where Mi(0) and M2(0) are the moments of the initial distribu- 
tion. Thus the variance a 2 = M 2 — Mf increases quadratically in time for 
arbitrary initial conditions. This result holds for arbitrary (non-trivial) values 
of the parameter 9. In Fig. 2 we show the time evolution of the relative differ- 
ence Aa 2 /a 2 where Act 2 = a 2 ap — a 2 is the difference between the approximate 
variance a 2 p obtained from eqs.(29) and (30) and the exact variance a 2 from 
the original Hadamard map, eq.(4). 

The independence of the quadratic increase of a 2 on the initial conditions has 
been established numerically for the particular case of the Hadamard walk [14]. 
Here we have also demonstrated the independence of this quadratic increase 
on the parameter of the unitary transformation. This is a new result. 



4 Conclusions 

The quadratic increase in time of the variance for the discrete time gener- 
alized quantum walk has been obtained analytically. This increase remains 
quadratic for arbitrary initial conditions and for all non-trivial values of the 
parameter 9 controlling the unitary transformation in the chirality subspace. 
This quadratic increase results from interference effects and thus it is strongly 
dependent on the coherence of the quantum evolution. If, for any reason, the 
quantum evolution is decoherent then the increase in the variance becomes 
linear with time. 

The general approach of separating the quantum evolution equation into a 
master equation supplemented by a term which takes into account quantum 
coherence effects provides a new perspective which is helpful in clarifying 
why the quantum evolution spreads faster than in the classical one. In the 
quantum case, there is a superposition of a left-propagating wave and a right- 
propagating wave, with both wavefronts traveling with constant speed cos 6*. 
Thus, their separation increases linearly in time and the variance does so 
quadratically. In the Markovian approximation, as developed in this work, the 
particle moves a step either to the right or to the left, making its choice in 
a random way. Due to the randomness of this motion, the variance increases 
only linearly with time. This process can be visualized as a frequent posi- 
tion measurement process, which amounts to reinitialize the system, at each 
measurement, in a different state. The wave function collapse causes memory 
loss of the previous distribution. Thus, we have a Markovian process which 
describes a series of random steps. 

We have shown how the Markovian approximation method can be applied 
to a generalized quantum walk, for arbitrary values of the parameter 9 and 
arbitrary initial conditions. When the evolution becomes decoherent, the in- 
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terference terms may be dropped and the evolution is described by a master 
equation implying a linear increase of the variance with time. The diffusion 
coefficient is not, in general, the same as that for the classical random walk. 
Only in the particular case of the Hadamard walk the diffusion coefficient is 
1/2 as in the classical case. This formalism shows in a transparent form that 
the primary effect of decoherence is to make the interference terms negligible 
in the evolution equation and then the Markovian behavior is immediately 
implied. This is true for any evolution operator of the form given in eq.(l). 

We have established the analogy between the generalized quantum walk and 
the resonant kicked rotor. This is done by obtaining an expression for the time 
evolution of the chirality amplitudes of the quantum walk and showing that 
they have the same form as the angular momentum components of the wave- 
function of the kicked rotor in a resonant regime. We related the kicked rotor 
strength parameter to the parameter 6 defining the unitary transformation of 
the chirality. This analogy opens several interesting possibilities. A quantum 
algorithm implements the evolution of the quantum kicked rotor in a quan- 
tum computer is known [18] . Since the quantum random walk on the line has 
the same long-time dynamics as the resonant kicked rotor, the same quantum 
algorithm may be used to describe the evolution of a generalized discrete time 
quantum walk on the line. The quantum kicked rotor has been experimentally 
realized using ultra-cold atom traps and some experiments have focused on 
the resonant case [19]. This opens interesting possibilities for the experimental 
realization of the quantum walk using quantum optics. 

Finally, note that the Markovian approach presented here provides a system- 
atic way to find the classical analog associated to a given quantum random 
walk. By reformulating the quantum problem as described here and neglect- 
ing the interference terms, the Markovian equation describing the equivalent 
classical problem can be readily obtained. 
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